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This Preliminary Amendment cancels without prejudice 



original claims 1-29 in the underlying PCT Application No. 
PCT/DE2004/002539 and adds new claims 30-58. The new claims 
conform to U.S. Patent and Trademark Office rules and do not 
add new matter to the application. 



Substitute Specification (including the Abstract, but 
without the claims) contains no new matter. The amendments 
reflected in the Substitute Specification (including 
Abstract) are to conform the Specification and Abstract to 
U.S. Patent and Trademark Office rules or to correct 
informalities. As required by 37 C.F.R. § 1.121(b) (3) (ii) 
and § 1.125(c), a Marked Up Version Of The Substitute 
Specification comparing the Specification of record and the 
Substitute Specification also accompanies this Preliminary 
Amendment. Approval and entry of the Substitute 
Specification (including Abstract) are respectfully 
requested. 



includes an International Search Report, dated April 28, 
2005. The Search Report includes a list of documents that 
were uncovered in the underlying PCT Application. 



present application is new, non-obvious, and useful. Prompt 
consideration and allowance of the application are 
respectfully requested . 
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DEVICE AND METHOD FOR ANALYZING A SAMPLE PLATE 
Background Inf ormat ion Field of the Invention 

The present invention relates to a device for the analysis of a 
sample plate according to the type more closely opacified in the 
generic part of Claim l containing material samples , as well as a 
5 method for analyzing the sample plate. 

Background Information 

It is of major importance, particularly in the field of material 
sciences, chemistry and pharmacy, to discover and develop 
optimized substances and materials with respect to 

10 rc3pcctivc desired cases of application. In this connection, sensor 
systems represent a special field of application, which are a key 
technology having a steadily growing number of applications, both 
in industry and in the private sector. For, acnooro Sensors are used, 
for example, in technical process monitoring, in the area of 

15 environmental protection, in the field of medicine and in the motor 
vehicle field. A considerable quantity of development work is 
currently being put especially into the development of faster and 
highly sensitive sensors having a low cross-sensitivity. 

As a rule, up until now, the developments were limited to 
20 optimization or modification of known materials. However, there 
is the problem that, for certain areas of application of the sensor 
systems, there is a great requirement for new materials which 
cannot be sufficiently covered by conventional methods which are 
distinguished by the production of individual sensors and a 
25 subsequent sequential characterization. 

In particular, in the development of new types of sensitive 
materials or material combinations, it may be expedient to use 
methods from the field of combinatorial chemistry or so-called high 
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throughput methods. In these methods, parallelized synthesis 
methods and screening methods are involved, by which new materials 
and material combinations may be discovered, or known synthesis 
methods for existing materials may be optimized in a broad field 
5 of parameters. 

A general representation of high throughput methods is known 
f ro m described in US Patent 5,925,355, and in thio which document— 
especially the application of combinatorial chemistry, known 
coocntially from in the field of pharmacy, to chemical and material 
10 science application areas bcing is proposed. 

A device of the type named in the — int roduc t ion , — for analyzing a 
sample plate is kno wn de scribed , for example, fro m in published 
German patent document DE 101 31 581 A3r. This device includes a 
sample plate on which 64 material samples have been applied in 
15 matrix fashion, which, in each case, are connected to two 

electrodes that, in turn, are provided with contact locations to 
which a means for reversible and addressable contacting to the 
installation may be brought. 

Summary of the Invent ion Summary 

2 0 The device— according to the present invention— for analyzing a 
sample plate having the fcaturco a ccording to the generic part of 
Claim l, — and having has a measuring head— that may be positioned 
in a housing carrier, which measuring head includes two measuring 
wires— per material sample for the electrical connection to the 

25 contact means, which lie against contact surfaces of the sample 
plate in a prestressed manner— and are connected to a measuring 
unit and an evaluation unit— - The device according to the present 
invention has the advantage of being simple to handle, since, in 
the use of a normalized sample plate, the contact between the 

30 material samples and the measuring unit and evaluation unit is able 
to be produced by simply setting the measuring head into the housing 
carrier. 
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The device according to the present invention is particularly- 
suitable for developing and discovering materials or material 
combinations that may be used as sensor materials, and may be 
characterized by their electrical properties. For example, the 
5 device according to the present invention may be used for 
developing an optimized sensor material of a gas sensor. 

Using the device according to the present invention, it is possible 
to almost simultaneously ^e-investigate a large number of potential 
sensor materials, disposed on the sample plate, under various test 
10 gases at different temperatures^ which may amount to up to 800°C. 
The investigation may be made by potentiometric , resistive, 
capacitative methods and also by complex impedance spectroscopy. 

The concept of spectroscopy should be understood, in this case, 
as involving frequency- dependent measurements, that is, the 

15 impedance of a sample is investigated at different measuring 
frequencies. For example, the individual material samples are 
tested in each case in a frequency range between 10 Hz and 10 7 Hz 
at a measuring data density of 15 measuring points per decade. This 
means that, per material sample, 180 measured data are ascertained 

20 at such a measuring data density. A far-reaching data reduction 
may be achieved, for instance, by the adaptation of a suitable 
circuit equivalent to the measured data. 

In order to hold the measuring wires which-r — in particular, lie 
against the contact surfaces of the sample plate via a fusion ball # 

2 5 with prestressing, the measuring wires may each be connected to 

an especially gilded spring contact, which will ensure a constant 
contact pressure of the respective measuring wire against the 
respective contact surface. 

In order to be able to apply a certain test gas to the material 

3 0 samples, for the development of a gas sensor, the measuring head 

may be connected to a gas supply unit. 
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In order. to be able to expose the material samples to different 
test gas atmospheres or reference gas atmospheres, the gas supply 
unit, which is expediently connected to a data processing unit of 
the measuring and evaluation unit, includes a gas mixing device. 
5 Furthermore, the gas supply unit may include a water reservoir for 
moistening the test gas atmosphere or reference gas atmosphere. 

Furthermore, the measuring head may be designed in such a way that 
it includes a gas chamber, as an integrated component, that is 
situated above the material samples of the sample plate and is 
10 preferably formed by an essentially bell -shaped distributing 
device. The gas chamber is connected to the gas supply unit. 

In order to achieve a homogeneous distribution of the test gas or 
the reference gas in the gas chamber, in one advantageous 



15 invention, a diffuser is disposed in the gas chamber. 

In the case of a plurality of material samples on the sample plate, 
in order to be able to measure the individual material samples in 
a simple manner, the measuring and evaluation unit advantageously 
includes two relay switch panels that are connected to the 

2 0 measuring wires, and, for example, in the case of 64 material 

samples on the sample plate, preferably each have a 3x64 matrix 
of relays appropriate for high frequency. In this case, the 64 
material samples on the sample plate are able to be measured in 
a measuring cycle, at least three measuring variables being 

2 5 accessible, namely, for example, the impedance via an impedance 
analyzer and direct current resistances of the material samples 
via additional appropriate measuring devices and their 
current /voltage characteristics curves . 

The measuring and evaluation unit is expediently equipped with a 
30 measuring and control software, which, on the one hand, controls 
the measuring progress and, on the other hand, passes along the 
measured data obtained to an appropriate data file or even to a 

NY01 1175944 vl 4 MARKED-UP VERSION OF 




ixample embodiment of the device according to the present 



SUBSTITUTE SPECIFICATION 



r 



relational databank, which may be read in by an evaluation 
software . 

The evaluation software preferably works in such a way that it 
includes a fit functionality for computing theoretical impedance 
5 spectra for the individual samples, the computation preferably 
being made based on a circuit equivalent that includes at least 
one virtual or real electronic component. A virtual component is, 
for instance, a constant phase element (CPE) . Thus, the fit 
functionality computes, based on the measured data for a circuit 

10 equivalent made up, for example, of a serial RC element, a 

theoretical impedance spectrum approximated as well as possible 
to a measured impedance spectrum, a variation of the capacitance 
and/or the resistance of the components of the RC element being 
carried out for the adaptation of the theoretical to the measured 

15 spectrum. Thus, in this case, a resistance value and a capacitance 
value for the RC element are assigned to the theoretical impedance 
spectrum that is adapted as well as possible. From the resistance 
value, for example, one may draw conclusions on the sensitivity 
of the respective material sample by comparing it to a reference 

2 0 value. 

And, — iftln the case of a plurality of material samples, which are 
measured under various measuring conditions, in order to make 
accessible to a simple evaluation the output values obtained for 
the individual material samples representing, for instance, the 
25 sensitivities, the evaluation software advantageously includes a 
data^mining functionality. The data^mining functionality 
determines the optimal material sample for the respective 
application case by a numerical path, for instance, by 
advantageously applying multidimensional object functions. 

30 Alternatively^ or in addition, the data^mining functionality may 
also include a visualization functionality. In this case, a user 
is able to determine the optimal material sample for the respective 
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application case, supported by a video screen. For instance, the 
visualization functionality works with a color spectrum, one 
certain color having a high sensitivity of the material sample 
assigned to it and another color having a lower sensitivity 
5 assigned to it . 

In order to be able also to discover sensor materials that are 
suitable for application cases in which, such as for instance in 
exhaust gas, where high temperatures prevail, the device according 
to the present invention preferably has a heating device into which 
10 the sample plate may preferably be immersed. In this case, the 
device especially has a**a high temperature reactor that is bordered 
by the heating device, and in which the sample plate that includes 
the material sample may have applied to it different test gases 
or reference gases. 

The present invention also te sprovides a method for analyzing a 
sample plate having the features of Claim 19 ao oubject matter . 
By applying this method, it is possible electrically to measure, 
under different conditions, a sample plate on which a large number 
of material samples, for instance, 64 material samples are disposed, 
and to select in a fully automatic fashion a material sample best 
suited to the application case. This takes place by an automatic 
selection of starting values for the components of the respective 
circuit equivalents and the subsequent error minimization 
computation. The starting values are applied in the error 
minimization computation, and starting from the starting values, 
a theoretical impedance spectrum, adapted to the respectively 
measured impedance spectrum, for the respective material sample 
being calculated under the respective measuring conditions. By 
this procedure, a considerable data reduction is made possible, 
since, starting from the impedance spectra which have a plurality 
of measuring points, few fit values representing derived variables 
are ascertained that describe the individual components of the 
circuit equivalents. The fit values thus represent dimensionings 
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of the components of the circuit equivalents, using which the 
measured impedance spectrum may be simulated as best as possible. 



In the determination of the circuit equivalents, for example, a 
serial connection of four RC elements may be selected, the higher 
RC elements, if necessary, being set to a starting value that has 
no influence on the computation of the theoretical impedance 
spectrum. 

The starting values of an RC element, required for the error 
minimization computation, are preferably computed from the 
maximally measured, imaginary impedance Z"_MAX and the 
corresponding measured frequency f_Z"__MAX, according to the 
following formulas : 

Rl START = -2-Z"_MAX 

1 



CI START = 



2x-f_Z"_ MAX • R _ START 



If the circuit equivalent includes several RC elements, the 
starting values for the error minimization computation required 
for the higher RC elements are preferably ascertained from the 
difference spectra between the measured data and data which are 
computed or simulated based on the starting values computed for 
the first RC element . 

In particular, the selection of so-called "good" starting values, 
which are close to the actual variables of the components of the 
circuit equivalent, decisively shortens the duration of the 
subsequent error minimization computation. "Bad" starting values, 
on the other hand, may lead to the error minimization computation, 
carried out based on the starting values, yielding meaningless 
values . 

In one preferred opocif ic example embodiment of the method 
according to the present invention, the circuit equivalent for the 
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simulation of impedance spectra, based on fit values obtained in 
response to an error minimization computation, is made up of a 
serial separation of four RC elements. By the determination of 
starting values for the individual components, it is then 
5 determined, in the light of a threshold value, how many RC elements 
are taken into consideration in the simulation computations. The 
threshold value is preferably a value that may be preset by a user. 
The percentage ratio of the resistance of the first RC element to 
the resistance of the current RC element n is checked according 
10 to the formula 

RCn_START > Wert [%] • RC1_START. 

The variable "value" is the variable, changeable by the user, which, 
for instance^ is preset to 10%. If the argument is not satisfied, 
the starting values for the components of the respective RC element 
15 are set to values that have no influence on a simulation of the 
impedance spectra. These values are kept constant in the error 
minimization computation. 

Furthermore, in the method, the validation variable is determined 
which evaluates the agreement between the calculated, theoretical 
20 impedance spectrum and the respectively assigned, measured 
impedance spectrum . 

The evaluation quantity represents the output value relevant to 
the respective analysis which, for example, in the determination 
of a sensor material for a gas sensor, reflects a sensitivity of 
25 the respective material sample. The sensitivity is a measure for 
the quality of a sensor. 

The sensitivity of a resistive gas sensor may be defined in various 
ways. If, for instance, one takes into consideration the direction 
of the change of the resistances during the application of gas, 
3 0 a sensitivity S may be expressed as the quotient of the resistance 

NY01 1175944 vl 8 MARKED-UP VERSION OF 

SUBSTITUTE SPECIFICATION 



R_TEST under a test gas atmosphere and the resistance R_0 under 
reference conditions, as follows 

R TEST . 

S = = ; for oxidizing gases 

R_0 

S = + — ^ -® — • f or reducing gases 
R _ TEST 

The sign of sensitivity S gives information with regard to the 
resistance change . 

Alternatively, the sensitivity may be described as a change in the 
resistances, namely, according to the formulas 

p A R TEST-R 0 , 

S _ A = = =— ; for oxidizing gases 

R _ TEST 

0 A R 0-R TEST ^ 

S _ A = + — = = ; for reducing gases. 

R 0 



These formulas yield sensitivities S_AA between -1 and 1, that is, 
normalized sensitivities . 

It is also possible to mutually— £e> convert- the sensitivities of 
the two definitions into each other. The sensitivity expressed by 
the value S is especially meaningful in the case of a large change 
of the resistance of a material sample that results from the 
application of a test gas. The sensitivity expressed by S_AA is 
particularly meaningful if the resistance of a material sample 
changes only a little as a result of the application of a test gas. 
However, the sensitivity expressed by S_AA has the advantage of 
a large tolerance with regard to measuring accuracies, a precise 
illumination of a region of small resistance changes, and therewith 
a better valuation of cross-sensitivities, as well as the 
possibility of an automated visualization and data processing. 
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In one particularly simple opecif io example embodiment of the 
method according to the present invention, the circuit equivalent 
is made up of a virtual arrangement of real electronic components 
such as capacitors and resistors. In its simplest embodiment, the 
5 circuit equivalent is composed of one resistor. 

By the error minimization computation, in which an adaptation of 
the variables of the components to the measured data takes place, 
a simulated impedance spectrum is determined for the circuit 
equivalent for the respective material sample. From the adapted 

10 variables of the components, one is able to draw conclusions on 
the electrical properties of the processes in the material sample 
described by the circuit equivalent. If the electrical behavior 
of the material sample is determined by several processes of 
different relaxation times, it is necessary to draw upon 

15 complicated circuit equivalents which, for instance, are made up 
of several RC elements connected in series. If the processes react 
differently to a variation of the measuring conditions, one may 
undertake an assignment of individual processes to components or 
groups of components so that the individual processes may be 

20 analyzed separately. Because of the error minimization computation, 
the data volume for the description of the impedance measurement 
is reduced to the arrangement and the magnitudes of the components. 

In the case of a sample plate having, for example, 64 material 
samples and a measurement under eleven different gas atmospheres 

25 at four different temperatures in each case, the quantity of 
measured data is approximately 2816 impedance spectra. In the 
method according to the present invention, the individual 
impedance spectra are simulated and in each case reduced to the 
derived measured variables that are reproduced by the variables 

30 of the components of the circuit equivalent. 

The computation of a theoretical impedance spectrum of a circuit 
equivalent is made byin one preferred opecif ic example embodiment 
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of the method according to the present invention in such a way that 
a complex admittance Y* and a phase shift cpcp of the individual 
components of the RC elements are determined at a given angular 
frequency coco (2Fin measuring frequency) according to the following 
5 formulas : 




0=arctan (coRC) 
Y'= cos0 
Y"= sin0 

10 The respective impedances may be assessed by a transformation: 

Y' Y" 

Z'= - : 

where Z' is the real part of the impedance and Z" is the imaginary 
part of the impedance . 

15 By a serial connection of the RC elements, the impedances of the 
individual RC elements may be summed directly. A calculation of 
the impedances for frequencies that correspond to the measured 
frequencies yields a data set corresponding to the measured data, 
so that an error estimation between the measured data set and the 

20 theoretically determined data set is possible. 

The error minimization computation carried out in the method 
according to the present invention is carried out in one 
advantageous embodiment by variation of the variables of the 
individual components by 1% . An error may be determined by analysis 
25 of the differences between the theoretically calculated spectra 
and the measured ones. If, after a variation, an error decreases, 
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a renewed variation is carried out for the same component of the 
circuit equivalent, and, based on this variation, a theoretical 
impedance spectrum is computed, using which a new error calculation 
is carried out. If the error does not decrease, the variation that 
5 has taken place is reversed, and the variation of the components 
having the opposite sign is varied, or another component is varied. 

The error computation is preferably carried out in such a way that 
the function for determining the error takes into account the 
starting value of the resistance of the first RC element. If the 

10 value is greater than a nominal measured resistance of the 

impedance analyzer used in the measurement, for instance, greater 
than 3- 10 7 QQ, the error is determined only from the imaginary 
part of the calculated impedance. In order to be able to weight 
the high-frequency range of the ascertained spectrum, the 

15 ascertained deviations are preferably multiplied by the logarithm 
of the measured frequency. This weighting makes it possible to 
suppress physically meaningless or even erroneous measured data 
in the low- frequency range. 

For the calculation of the validation variable, that is, for the 
20 validation of the quality of the theoretically calculated 
impedance spectrum for the respective material sample, 
prcf crably e .g. , a corridor is determined about the theoretical 
impedance spectrum which contains, for example, 90% of the measured 
data. In order to minimize the time required for the computation 
25 of the validation variable, a successive approximation algorithm 
may be used for this, the boundary values of the algorithm 0 QQ 
and of the double values of the real parts of the summed impedance 
may be at the smallest measured frequency. 

The method according to the present invention io prof Grably m ay be 
30 used for the analysis of material samples under various test gas 
atmospheres and especially at various temperatures. In this case, 
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impedance spectra are measured under various measuring conditions 
for all material samples situated on the sample plate. 



In the method according to the present invention, the respectively 
measured spectra are then each simulated in a theoretical manner 
5 by an error minimization computation based on a circuit equivalent . 
This results in a large quantity of evaluation magnitudes which 
represent the target variables in the method. 

In the method according to the present invention, the target 
variables or the evaluation variables are preferably written into 
10 a databank and evaluated using a datamining functionality. 

As far as the databank is concerned, which is used for receiving 
and making available data sets, expediently a relational databank 
is involved in which the data are filed in the form of tables, 
ordered according to thematic areas. Relationships between the 
15 individual data sets in the tables are produced by so-called 
identification keys. 

The databank may, for example, include additional properties of 
the material samples, such as the synthesis conditions of their 
starting materials, their sampling history and the like. These 
20 properties are linked with one another via table relationships. 

The datajMTiining may be performed using a possibly 
multidimensional target function and/or using a visual 
data^mining functionality. Data^ming, carried out using a target 
function, is a numerical method which is based on the individual 

25 evaluation variables which are stored in the databank, for instance, 
classified according to measuring temperature and test gas. In this 
context, it is first specified which properties are wanted of the 
material being sought. For instance, in the search for a sensor 
material for a gas sensor it may be stated with respect to which 

30 gas the sensor material is to be sensitive, and which 

cross-sensitivities could have an interfering effect. Thus, from 
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this there develops a demand profile for a fingerprint of the 
sensitivities . 

The visual data_^mining functionality advantageously works in such 
a way that the evaluation variables of the material samples of the 
5 sample plate are represented broken down, for instance, according 
to test gases and temperatures. 

In the use of the device according to the present invention and 
the method according to the present invention, one has a preferably 
fully automatic high- throughput impedance system using which new 

10 materials having a high sample throughput and a low time 

expenditure and cost expenditure are able to be developed. Thus, 
for example, by using the system on two days, 64 different material 
samples may be investigated and also evaluated at four different 
temperatures and under eleven different test gas atmospheres, with 

15 regard to their sensoric properties. 

In particular, the system may be used for general material 
development and especially for the development of sensors in the 
motor vehicle field and in the safety technology field. 

The present invention also h as encompasses data processing 
20 equipment having a data processing program as subject matter, for 
implementing the method according to the present invention. 

Further advantagco and advantagcouo rcfincmcnta of the subject 
matter in accordance with the prcocnt — invention rcoult — from the 
description, — the drawing and the patent claims. 

25 Brief Description of the Drawings 

An exemplary embodiment of the device according to the present 
invention io shown in a oimplificd version in the drawing and io 
elucidated in more detail in the following description, — if* 
connection with a method according to the present invention. . — The 
3 0 figures show: Brief Description of the Drawings 
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Figure 1 shows a test setup of a device according to the present 
invention^ 

Figure 2 shows a measuring appliance of the device according to 
Figure dhrl . 

5 Figure 3 shows a measuring head of the measuring device according 
to Figure 2-^2 . 

Figure 4 shows a gas distribution device of the measuring head as 
in Figure 3-^3 . 

Figure 5 is a flowchart illustrating a measuring sequence in the 
10 light of a flow diagram; 

Figure 6 shows a diagram representing a measuring sequence over 
time for material samples at various temperatures and under various 
gas atmospheres-t-^ 

Figure 7 shows a diagram in which an influent behavior is shown 
15 e #for four different surface doped material samples^ 

Figure 8 shows a diagram in which a response time performance is 
shown e #for three different material samples during a test gas 
pulse-t-^ 

Figure 9 shows a circuit equivalent^ 

20 Figure 10 shows a flow diagram of an error minimization 
c ompu t a t i on^ 

Figure 11 shows an exemplary curve of the sensitivity of a material 
sample^ 

Figure 12 is a graph illustrating an exemplary fingerprint of 
25 sensitivities of a material sample with regard to various gas 
atmospheres-? — and_ ._ 
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Figure 13 shows a greatly schematized representation of a visual 
data^mining functionality . 

Detailed Dcocription of the Exemplary Embodimcnt Detailed 
Description 

5 In Figures 1 through 4, a test setup is shown of a device 10, for 
analyzing a sample plate 12, on which 64 material samples 13 are 
situated. Device 10 includes a high temperature reactor 14, a gas 
supply unit 16 and a measuring and evaluation unit 18. 

High temperature reactor 14, which is shown in more detail 
10 especially in Figures 2 through 4, includes a frame or housing 20, 
in which an adjustable-height heating block 22 is situated, which 
is supported on a guide rod 23 and a threaded rod 24 that is provided 
with a crank 25. 

Heating block 22, which is made up, for instance, of four heating 
15 plates that border on a heating chamber, having a heating power 
of 1100 W each, is adjustable in height in such a way that plate 
carrier 28, and thus sample plate 12^ may be immersed into the 
heating chamber. 

Furthermore, high temperature reactor 14 includes a measuring head 
20 26 which is set into measuring head carrier 27 that is connected 
to frame 20, and is used for contacting the 64 material samples 
13 situated on sample plate 12 to measuring and evaluation unit 
18. Sample plate 12 is situated on a plate carrier 28 that is 
connected to measuring head carrier 27. 

25 Measuring head 26, which is particularly shown in Figure 3, 

includes a ground plate 29, which is formed, for example, from a 
workable glass ceramic and is used as a mounting support of 12 8 
measuring wires 30A, 30B, that are each made of platinum and are 
sheathed by an aluminum oxide tube. At their lower ends, measuring 

30 wires 30A, 30B each have a fusion ball 31A, 31B which are used for 
contacting each measuring wire to a contact surface of sample plate 
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12. Sample plate 12 has two contact surfaces per material sample, 
i.e., in the present case altogether 128 contact surfaces, which 
each collaborate with one of measuring wires 3 OA, 3 OB. Measuring 
wires 3 OA, 3 OB are provided at their upper ends, in the region of 
5 ground plate 29, each with a gilded spring contact 32A, 32B, which, 
when measuring head 26 is in a fixed position, ensures a constant 
contact pressure of measuring wires 30A, 30B and of fusion ball 
31A, 3 IB against respectively assigned contact surfaces of sample 
plate 12. Furthermore, the sheathings of measuring wires 30A, 30B, 
10 that which are made of platinum, are fixed to supporting plates 33 
and 34, which in turn are fastened to a rod assembly 35, and are 
aligned parallel to ground plate 29. 

Spring contacts 32A, 32B are connected via lines 36A, 36B to SMB 
sockets 37A, 37B that are integrated into measuring head 26, and, 

15 in turn, are connected via shielded SMB lines 3 8A, 3 8B to measuring 
and evaluation unit 18 of device 10. Altogether, measuring head 
26 has 128 SMB sockets 37A, 37B, of which each is connected to a 
measuring wire 30A, 30B and to which, in each case, one SMB line 
38A, 38B is connected, that leads to measuring an evaluation unit 

20 18. However, for clarity's sake, in Figures 1 to 3, in each case, 
only two SMB sockets 37A, 37B and two SMB lines 38A, 38B are shown, 
as well as measuring wires 30A, 30B that are allocated in each case. 

Furthermore, measuring head 26 has a gas distribution device 39, 
made, for example, of quartz glass and essentially bell-shaped, 
25 which is shown in detail in Figure 4 and which is connected via 
a gas supply line 40, made of high-grade steel, to gas supply unit 
-1-8-16 . Bell -shaped gas distribution device 3 9 borders on a gas 
chamber that is situated above the 64 material samples 13 that are 
disposed on sample plate 12 . 

30 Material samples 13 are formed, for instance, of tin oxide Sn0 2 , 
and have different dopings which are formed, for example, of 
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lanthanides. Material samples 13 are distributed in a matrix-like 
manner, in eight rows and eight columns, on sample plate 12. 

Gas distribution device 39, that is shown in detail in Figure 4, 
also has in its gas chamber 41 a diffuser insert 42, formed of a 
5 quartz ball and has a plurality of bore holes 43 that each have 
a diameter of about 1 mm. 

In order to manage a uniform gas outflow from gas chamber 41, 
bell-shaped gas distribution device 39 has at its edges spacers 
43- 4 3 0 , which establish a gap of 0.8 mm in width between gas 
10 distribution device 39 and sample plate 12. 

In order for gas chamber 41, and therewith material samples 13, 
to have different test gas atmospheres applied to them, gas supply- 
unit 16 has two gas bottle cabinets, not shown here in detail, each 
having four gas bottles, which are respectively connected to a gas 

15 flow controller 44A, 44B, 44C, 44D, 44E, 44F, 44G and 44H, one gas 
bottle containing moist synthetic air, the second gas bottle 
containing hydrogen, the third gas bottle containing methane, the 
fourth gas bottle containing synthetic air, the fifth gas bottle 
containing nitrogen dioxide, the sixth gas bottle containing 

2 0 nitrogen monoxide, the seventh gas bottle containing 

propanc propene and the eighth gas bottle containing carbon 
monoxide. The capacities of flow controllers 44A, 44B, 44C, 44D, 
44E and 44F are respectively between 0 seem and 100 seem. The 
capacities of flow controllers 44G and 44H are respectively between 

2 5 0 seem and 10 seem. By appropriate control of the. volume flow of 

the various gases, using the group of uniform devices, coming from 
the eight gas flow controllers 44A, 44B, 44C, 44D, 44E, 44F, 44G 
and 44H test gases of different compositions may be fed into a 
collecting line 45 that is connected to gas supply line 40. 

3 0 In order to set a relative humidity of the respective test gases, 

the test gas may have admixed to it a moist carrier gas, consisting, 
for example, of synthetic air. The moisture in the carrier gas is 
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set by bubbling it through a water reservoir 46. Measuring the 
humidified carrier gas using a moisture sensor, not shown here in 
greater detail, yields, for example, a relative humidity of 
approximately 90% at room temperature. 

5 Measuring and evaluation unit 18 includes two relay switch panels 
50 and 51, to which respectively 64 lines 38A, 38B are connected 
that lead to measuring head 2 6 of high temperature reactor 14 . Relay 
switch panels 50 and 51 each form a 3x64 matrix of relays suitable 
for high-frequency. 

10 Relay switch panels 50 and 51 are connected to a measuring and 
evaluation computer 53 via a digital control line 52. Moreover, 
relay switch panels 50 and 51 are connected to an impedance analyzer 
64 and a so-called source meter 55, via measuring lines 54. These 
two measuring unite are also connected to measuring and evaluation 

15 computer 53 via digital control line 52. The addressing of 

impedance analyzer 54 and source meter 55 takes place via the two 
relay switch panels 50 and 51. 

Measuring and evaluation computer 53 is also connected via an 
additional digital control line 56 to a D/A-A/D converter 57, which 
20 is connected via an analog control line 58, on the one hand, to 
heating dcvicc block 22 of high temperature reactor 14, and, on the 
other hand, to gas flow controllers 44A, 44B, 44C, 44D, 44E, 44F, 
44G and 44H of the gas supply unit _[or gas mixing battery_)_ 16. 

On measuring and evaluation computer 53, a modular measuring and 
25 control software is filed, which makes possible, via a script 

control, a complete automation of measurements carried out using 
device 10. Because of the modular design of the measuring software 
and control software, a broadening of the measurement system is 
possible without a problem. 

30 Using source meter 55, direct current resistances, U/I 

characteristics lines or voltages of individual material samples 



NY01 1175944 vl 



19 MARKED -UP VERSION OF 

SUBSTITUTE SPECIFICATION 



13 may be measured. These values, as well as the measured values 
ascertained using impedance analyzer 54, may be passed on via 
digital control line 52 to measuring and evaluation computer 53, 
which includes a databank for the measured data. 

5 In order to measure material samples 13, that are disposed on sample 
plate 12, in a high throughput mode or according to a high 
throughput method, via a script data file is used, a list of tasks 
that is made up of index key words and parameters being passed over, 
via a script data file, to the software^ using which the script 

10 data file is processed and which controls all the functions of the 
system. The script data files ensure a continual checking of 
control parameters, so that measurements under faulty measurement 
conditions are excluded in that a further processing of the script 
data file is suspended. There is no limit with regard to the length 

15 of the script data file. 

A measuring sequence able to be carried out using the device 
according to Figure 1 is shown in Figure 5 in the light form of a 
flow chart. For the analysis of resistive gas sensor properties, 
material samples 13 of sample plate 12 are first electrically 

2 0 characterized under a reference gas atmosphere which, for example, 

is formed of synthetic air of a relative humidity such as 45%, and 
under various test gas atmospheres. By a modulation of the 
measuring temperature, data may be gathered on the influence of 
the respective operational temperature on the material samples 
25 that are suitable as sensor material, as well as on the activating 
energy of conductivity processes. 

In order to ensure comparability of measured data of several sample 
plates, the script data file used is designed as standard script 
which controls the complete high throughput screening, whose 

3 0 chronological sequence is described in the flow chart shown in 

Figure 5. In this instance, in a first method step Ml, a measuring 
temperature T_MESS of sample plate 12 is set. Measurements may be 
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carried out, for example, in a temperature range between 400° C and 
250° C, in steps of 50°. During the cooling-off phases, sample plate 
12 has applied to it a reference gas of synthetic air having a 
relative humidity of 45%. 

5 When there is a change in the sample temperature, a stable base 
resistance of the material samples first sets in after a certain 
time, so that conditioning of the material samples in a step M2 
is required. Temperature changes in the material samples lead to 
a metastable state of intrinsic defects, which may represent oxygen 

10 vacancies, and whose thermodynamic equilibrium setting may require 
a finite time. Figure 7 shows in an exemplary manner the 
conditioning behavior or influent behavior for four different 
surface-doped ln 2 0 3 samples upon reaching a target temperature of 
300° C. Resistances R of the samples approach a boundary value, 

15 asymptotically with time, which represents the so-called base 

resistance or reference resistance. This boundary value is reached, 
as extrapolated, after about 90 minutes. To ensure a constant 
reference resistance, it is of advantage to select a conditioning 
time of 120 minutes. 

2 0 Subsequently, the test gas required for the measurement is put 

together in a method step M3 , and is introduced into gas chamber 
41 above sample plate 12 . The first test gas includes, for instance, 
hydrogen at a concentration of 25 ppm. Synthetic air acts as the 
carrier gas. The gas flow is set to 100 seem using moist synthetic 
25 air. 

In order to ensure that the individual material samples have 
reached their base resistance independent of their relative 
position on the sample plate, there subsequently takes place, in 
a step M4 , an advance phase for the test gas that was analyzed before . 

3 0 Figure 8 shows the resistance pattern of three material samples 

lying on a diagonal of the sample plate, during a test gas pulse 
lasting 4 0 minutes, the pulse gas containing propane at a 
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concentration of 50 ppm. The base material of the material samples 
is tin oxide Sn0 2 in the example shown in Figure 8, because of its 
great sensitivity to hydrocarbons. Independently of the relative 
position on the sample plate, the resistance of the material 
5 samples decreases, because of the test gas pulse, to a constant 
value, within approximately 6 minutes. Thus, the response time 
amounts to 6 minutes in each case, the response behavior being 
essentially independent of the position of the respective material 
sample on the sample plate. After terminating the test gas pulse, 

10 the resistances tend asymptotically to the base resistances within 
approximately 10 minutes. Accordingly, the response behavior is 
essentially independent of the position of the material sample on 
the sample plate. The advance time, selected for the script data 
file, of each test gas or reference gas amounts to about 15 minutes 

15 before the corresponding measurement. 

In a subsequent method step M5 , there takes place the measurement 
of impedance spectra for the 64 material samples disposed on the 
sample plate. In this context, in the script data file, for the 
measurements of the impedance spectra, the following are 
20 established as parameters: 



amplitude of the measuring voltage [V] 0.1 

starting frequency [Hz] 10 

ending frequency [Hz] 10 7 

measuring points per frequency decade 15 

25 bias [V] 0 

mode [HS: highs peed; NO: normal; AV: average] HS 



The measuring data obtained are ascertained using impedance 
analyzer 64 and passed on to measuring and evaluation computer 53 
or the databank filed on it . 

3 0 Thereupon the measurement may be carried out under another test 
gas atmosphere, in this instance, once again the program runs 
through steps M4 and M5 . For example, as further test gases, carbon 
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monoxide at a concentration of 50 ppm, nitrogen monoxide at a 
concentration of 5 ppm, nitrogen dioxide at a concentration of 5 
ppm or propane at a concentration of 2 5 ppm may be used. 

The entire screening may then be carried out by j umping looping back 
5 to step Ml at another measuring temperature . 

Figure 6 vioualizcG is a diagram that illustrates a measuring 
sequence determined by a standard script, in thio diagram, 
measurements under reference conditions and during conditioning 
phases not being shown— for the sake of clarit y in this diagram , 

10 Measurements under reference conditions, i.e.^ under a reference 
gas atmosphere, are always carried out before the initiation of 
a new test gas H 2/ CO, NO, N0 2/ or propanc propene at the appropriate 
temperature T. There is always respectively a so-called test gas 
advance phase X and a measuring phase Y. Furthermore, the 

15 respective concentrations C of the test gases are shown in Figure 
6. 

The measuring and evaluation software may have a functionality on 
which the respectively measured impedance spectrum is graphically 
shown on a monitor. The measured data, which arc profcrabl ^ ay be 
2 0 stored in ASCII format in the databank or in the index assigned 
to the sample plates on the measuring and evaluation computer, are 
consequently directly accessible to visual control. 

Furthermore, the measuring and evaluation software may be 
furnished with a functionality which shows the impedance spectra 

25 of the sample materials of a sample plate as an image matrix 

independently of position. In this context, raw data as well as 
derived data may be shown. The functionality may also include an 
evaluation window to which the data of a certain material sample 
may be transmitted. The functionality may also be combined with 

30 additional image -producing measuring or evaluation systems. Image 
data from these systems may then also be read into the image matrix 
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independent of position, in order thus to be able to obtain further 
data on the samples. 

After measuring the impedance spectra for the individual material 
samples, a theoretical impedance spectrum based on circuit 
5 equivalents is computed for each material sample. In the present 
case, the circuit equivalent is made up of a serial connection of 
four RC elements of the kind shown in Figure 9. The theoretically 
calculated impedance spectra, in this context, are each adapted 
as well as possible to the correspondingly measured impedance 
10 spectrum, namely, by varying the magnitudes or dimensionings of 
the individual components of the circuit equivalents. 

This is clarified in the following with reference to the error 
minimization computation in the light of the flow chart shown in 
Figure 10. 

15 In order to carry out the error minimization computation, it is 
necessary to determine starting values for the individual 
components of the circuit equivalents. At the same time, it is 
necessary to determine the number of RC elements to be taken into 
consideration in the error minimization computation. For these 

20 purposes, in a step SI, the maximum imaginary impedance Z"_MAX is 
ascertained from the measured data for the respective material 
sample. In the selection, it is also judged whether, in the case 
of the maximum measured imaginary impedance, a fault measurement 
is involved, which is done by an investigation of the imaginary 

25 impedances measured at adjacent measuring points, i.e.^ by an 
investigation of the local maximum. 

Subsequently, in a step S2 , starting values R1__START, for the 
resistance, and C1_START, for the capacity of the first RC element, 
are calculated based on Z"_MAX and the corresponding measuring 
3 0 frequency, according to formulas 

Rl START = -2 • Z" MAX 
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CI START = 1 



2n • f _z" _A4AX • r _start 

where f_Z u _MAX is the measuring frequency for the maximally 
measured, imaginary impedance. 

Subsequently, in a loop step S3 , starting values for the resistances 
5 and capacitances of the additional three RC elements of the circuit 
equivalents are ascertained. 

To do this, first of all, in a step S4 , based on starting values 
R1_START and C1_START, a theoretical impedance spectrum is 
calculated, and, in a step S5, a difference spectrum between the 

10 theoretical impedance spectrum and the measured impedance spectrum 
is calculated, and from this difference spectrum, in turn, the 
maximum of the imaginary impedance Z"_MAX is ascertained. Starting 
from this maximum, in a step S6, the starting values Rn_START and 
Cn_START (n = 2 to 4) are then calculated for the respective RC 

15 element being examined, according to the formulas explained in 
connection with step S2 . 

Thereafter, in a step S7, a threshold value examination is carried 
out for calculated resistance starting value Rn_START, using which, 
it is determined whether the respective RC element should be taken 
20 into consideration in the subsequent simulation calculations. The 
threshold value is a value that may be changed by the user. In the 
comparison, the percentage ratio of the resistance of the first 
RC element to the resistance of the current component is checked. 

If the calculated starting value Rn_START is greater than the 
25 threshold value, the RC element is taken into consideration, the 
system goes back to step S3, and a calculation of the starting 
values Rn_START and Cn_START is undertaken for the next RC element . 
If the argument is not satisfied, the magnitudes of the components 
of the higher RC elements are set to values that have no influence 
30 on a simulation calculation for an impedance spectrum. In 
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subsequent adaptation steps, these values remain constant, and 
they amount, for example, to Rn_START = 1 and Cn_START = 10~ 15 . This 
determination takes place in a step S8 . 

Thus, the threshold value query determines the number of the RC 
5 elements that are to be taken into consideration in the following 
steps. If the threshold value is undershot, the system transfers 
directly to the error minimization computation based on the number 
of RC elements determined, for instance, based on an RC element, 
in this case, a loop loop starting with step S9 being executed over 
10 all m RC elements taken into consideration. 

Within loop starting with step S9, an additional loop starting with 
step S10 is included , in which runs run 1 to 3 isare executed, in 
this case an empirical number of run-throughs being involved, which 
are carried out to increase the accuracy of the error minimization 
15 computation . 

Within loop S10, first of all, in a step Sll, resistance value 
R1_START of the first RC element is drawn upon as variable VAR. 
The variation of the individual components in the error 
minimization computation amounts to 1%, which is expressed by a 
20 value FAK = 0.01 in step Sll. The error is preset to 10 99 . The 
variation of the respective component takes place in a step S12, 
according to formula VAR = V--FAK. A calculation of a theoretical 
impedance spectrum then takes place, based on the varied value, 
in a step S13 . 

25 After that, an error is calculated which is based on a comparison 
between the impedance spectrum calculated for the correspondingly 
adapted component variables and the measured impedance spectrum. 
A function for determining the error is dependent on the starting 
value of the resistance of the first RC element. In order to 

30 calculate the error, in a step S14 it is checked whether the 

respective starting value is greater than the nominal measuring 
resistance of the impedance analyzer used, for example, greater 
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than 3 x 10 7 QQ. If this is so, the error is determined in a step 
SI 5, only by examining the imaginary parts of the impedance, 
according to the formula 

error =\ Z" _fit - Z"_mess | •log/ 

5 for the weighting of the high-frequency range of the spectrum a 
multiplication being made by the logarithm of the measuring 
frequency, whereby physically meaningless or erroneous measured 
data in the low- frequency range are suppressed. Otherwise, in a 
step S16 an error examination is made, based on both the real parts 
10 and the imaginary parts of the impedances according to the formula 

error = | Z" _ fit - Z" _mess | + | Z % _fit - Z' _mess | . 

In a step 17, the system checks whether the error has become smaller 
because of the variation of the component in step S12 . If this is 
the case, the respective component, for instance, the resistance 

15 of the respective RC element, is varied once more by going back 
to step S12 . If the error has not gone down, the variation that 
has taken place in step S12 is reversed in a step S18, and, in a 
step S19 it is established whether the sign of the variation will 
be changed in a step S2 0, or whether, in a step S21 and a subsequent 

20 step 22, the next component, in this case the capacitance C, will 
be selected as the variable component. In this case, then, there 
is once more a return to step S12, the magnitude of the component 
is then again varied as long as the error is decreasing. 

The variation of the resistance and the variation of the 
25 capacitance are repeated by return to step S10 twice for each RC 
element . Thereafter, as a function of the number of the RC elements, 
that are to be taken into consideration, there is an n-fold return 
to step S9. 

If no further minimization of the error is possible by variation 
3 0 of the resistance and the capacitance of the RC element, the values 
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for the resistance and the capacitance yielding the smallest error 
are output as evaluation variables of the error minimization 
computation . 

At the close of the error minimization computation taking place 
5 in steps S9 to S22 , in a step 23- S23 , a validation function is 
ascertained for the theoretically determined impedance spectrum, 
by which the quality of the theoretically calculated spectrum may 
be estimated. Upon determination of the validation function, for 
the estimation of a quality of the error minimization computation 

10 a corridor is ascertained, around the calculated impedance 

spectrum, that includes 90% of the measured data. The error thus 
ascertained lies between 0 and 1, an error of 0 expressing an ideal 
agreement of the measurement and the simulation calculation, and 
an error of 1 expressing no agreement of the measurement and the 

15 simulation calculation. 

When using the method for establishing a sensor material of a gas 
sensor, sensitivities S_AA are calculated, based on resistances 
of the circuit equivalents obtained by the error minimization 
computation. Sensitivities S__AA thus obtained are normalized, and 
2 0 lie between -1 and +1, as may be seen in the diagram in Figure 11. 
Figure 11 shows sensitivities for an assumed reference resistance 
of 100 £2Q, at varying test resistance R_test . Sensitivities S_AA 
are stored in the databank. 

Subsequently, based on the ascertained sensitivities S_AA, a 

2 5 numerical data^mining is carried out, so as to be able to select 

an optimal gas sensor for a special application case. The 
data^mining is carried out especially in a numerical way, the 
individual sensitivities S__AA being accessed which are stored in 
the databank, classified by measuring temperature and test gas. 

3 0 It is first specified which properties are demanded of the desired 

sensor material. In the simplest case, it is stated for which test 
gas the sensor is to be used, and which cross-sensitivities will 
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interfere. This yields a requirement profile on the fingerprint 
of the sensitivities, as shown in Figure 12. The requirement 
profile is converted into so-called or ^<^_ requirements, 

whereby data sets having the desired properties are characterized. 
5 Additional conditions, such as the type of base material of the 
sensor material or its surface doping may be taken up into the 
target function. In the example in Figure 12, the requirement 
profile turns out so that all test gases A, B, C, E, except for 
test gas D, have a sensitivity less than 0.2 and greater than -0.2, 
10 and test gas D has a sensitivity greater than 0.9. In each case 
X states the reference. Using an evaluation functionality, the 
requirement profile is converted to an SQL filter scan and passed 
on directly to the databank. The results of the SQL instruction 
may then be shown in tabular form on a screen. 

15 Alternatively, or in o upp lcmcntati on add i t i on , the sensitivities 
ascertained for the various material samples at the various 
measuring conditions may be evaluated using a visual data^mining 
functionality. In this functionality, which is shown in Figure 13, 
sensitivities of so-called library plates are shown classified as 

20 to test gases and temperatures on a screen. In the exemplary 

representation as in Figure 13, library plates 101 to 112 are shown 
for four temperatures A, B, C and D, and for three different test 
gases I, II and III. Each library plate 101 to 112 is assigned to 
one temperature A, B, C or D, and to one test gas I, II or III. 

25 Sensitivities S_AA of the individual materials of library plates 
101 to 112 are each shown according to their position on the sample 
plate as circles 120 and in an off -shade. In Figure 13, for the 
sake of clarity, only four of the 64 material samples are shown 
for each library plate 101 to 112. Positive sensitivities S_AA, 

3 0 for example, are shown in colors that go from black to yellow via 
red shades, whereas negative sensitivities are shown, for example, 
in colors that go from black to turquoise via blue tones. The 
diameter of individual circles 120 is determined by the validation 
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of the error minimization computation, in the case of a large error, 
a small circle being shown, and in the case of a small error, a 
comparatively large circle being shown. Using this scheme makes 
an intuitive judgment of large data quantities possible for the 
5 user. 

The data sets required for the functionality shown in Figure 13 
are extracted directly from the databank, so that the measurements 
to be shown may be selected via filter functions. The number of 
the data sets indicated is only limited by the working memory of 
10 the evaluation computer used. 

Furthermore, it is possible in the visual data^mining 
functionality that additional data, concerning the individual 
material samples, may be indicated by the screen- supported 
selection of a certain circle. 
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Abotract 



ABSTRACT 



A device and a method are provided for the analysis of a sample 
plate on which at least two material samples are situated. In the 
method, one impedance spectrum is measured for each of the material 
samples. As a function of the respectively measured impedance 
spectrum, a dcoign conf iguration of a circuit equivalent is 
determined which includes at least one electronic component. Then, 
for an error minimization computation, starting values for the 
components of the respective circuit equivalents are determined. 
In the error minimization computation, a theoretical impedance 
spectrum is calculated for at least one of the material samples, 
based on the impedance spectrum measured for the material sample, 
as well as the starting values for the components of the respective 
circuit equivalent, and fit values are determined for the 
components of the respective circuit equivalent. Subsequently, a 
validation variable is determined for the calculated, theoretical 
impedance spectrum, and an evaluation variable is ascertained by 
comparison of at least one of the fit values for the components 
to a reference value. 
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